NASA Tccfankai Memocandiim SM 


IMPflOVEMENTS TO THE lANGlEY HZE 
ION MODEL 



laivrence W. Townsend and HotI BL Bidasaria 


September 1982 



(M4Si*Ta>84542) IHPBOVEfiEMTS TO THE L4HGLET Ii6i- 13971 

HZE 43BAS10II BODEL (NASA) 21 p 

HC A02/HF AO I CSCL 20H 

Unclas 
G3/73 02045 




Nahonal Aeronautics and 
Space Admneiration 


Hampton. Virginia 23S6S 


Table or Tontents 


r^NCLATURE i 

SUHWRY 1 

INTROOUCTIOH 2 

ANALYSIS 3 

RESULTS 5 

CONCLUDING REMARKS 7 

REFERENCES 8 

TABLES 10 

FIGURES 11 


1 





NonencUturc 


A nuclear nass number 

a oscillator parameter, fm 

B(e) average slope parameter of nucleon-nucleon scattering amplitude, 
fm2 

t projectile Impact parameter vector, fm 

C average correlation function 

e two-nucleon kinetic energy In their center of mass frame, GeV 

kp Fermi momentum wavenumber, fm*’^ 

I(S) defined In equation (3) 

n number of abraded nucleons 

N neutron number 

position vector, fm 

r^ nucleon effective root-mean-square radius, fm 
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neutron root-metn-squere chtrge radius* fie 


proton rootHwan-square charge radius* ha 
two-nucleon relative position vector* fm 

total nueber of nuclear protons 

position vector of projectile In beaa direction* fm 

binomial coefficient 

collection of constituent relative coordinates for target* fa 
nuclear density* fm"* 

average nucleon-nucleon total cross section* fa^ or mb 
heavy-ion absorption cross section* fa* or ab 
cross section for abrading n nucleons* fia* or mb 
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Ii|>rovenents to the Langley HZE Abrasion Nodel 


by 

Laurence W. Toimsend, Langley Research Center, Kaaipton, Virginia 

and 

Harl B. Bidasarla, Old Dominion University, Norfolk, Virginia 


Summary 

Improvcnents to a previously developed HZE abrasion model are made by 
Incorporating more realistic values for the constituent Fermi momentun and 
nucleon root-mean-square charge radius. The theoretical predictions for 
neon projectiles at 2.1 GeV/nucleon colliding with carbon and molybdenum 
targets are In excellent agreement with recent experiment results. 



INTRODUCTION 


The attractiveness of HZE attenuatloi by nuclear fragmentation, as a 
means of radiation protection for future manned space aopllcatlons, 
dictates that a quantitatively accurate nuclear fragmentation model be 
developed. In previous work (refs. 1 and 2) an HZE abrasion model, which 
incorporates Pauli correlation effects and realistic density distributions, 
has been developed- In reference 2, the Importance of Pauli effects and 
the proper choice for the nuclear density distribution were clearly 
demonstrated by comparison with recent experimental results (ref. 3). For 
simplicity, the constituent Fermi momentuR chosen for use In that work was 
the value for Infinite nuclear matter. In addition, the nuclear 
distributions were obtained by unfolding the finite proton charge 
distribution from the experimental nuclear charge densities. In this work. 
Improvements In these two areas are made by utilizing more realistic values 
for the Fermi momentum (refs. 4 and 5) and by utilizing an "effective" 
nucleon charge distribution which accounts for the differences betweei the 
neutron and proton charge distributions (ref. 6). 
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ANALYSIS 


From reference 2, the cross section for abrading projectile nucleons 
Is 


>n *(!)*) / a - exp (-Ay o(e) 1(6)])" exp[-AyAp o(e) 1(6)] b db 


( 1 ) 


where the residual fragment (prefragment) mass number Is 

Ap - Ap - n (2) 

and 1(6) Is 

1(6) * [2 t » B(e)]"^^^ / dz / d^ly Pj(^y) I d^y Pp(6+z+y+|j) 

r ♦ , -v^ (3) 

[1 - C(y)l exp^gj 

The Pauli correlation function approximation C(y), Is 

C(y) * j exp (-k^ y^/10) (4) 

where kp » 1.36 fm"^ for Infinite -uc'ear matter. Since infinite nuclear 
matter Is approached only fcr very heavy nuclei (ref. 3), the corresponding 
Fermi momentum, although a reasonable approximation. Is generally an 
overestimate, especially for 
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lighter nuclei. This can be seen from Table I which lists values for Fermi 
momenta (kp) as a function of mass number, obtained from 500 NeV electron 
scattering experiments (ref. 4). The value for kp Is even smaller If the 
Incident lab momentum per nucleon (drift momentum) Is accounted for 
(ref. 5). 

The nuclear densities, pj and pp, shown In equation (3) are 
obtained by unfolding the gauss^sn nucleon charge density from the 
experimental nuclear charge distribution using the methods of references 2 
and 7. In those works, the nucleon charge density was assumed to be 
Identical to that of a bare proton. Since the charge distribution of a 
proton differs from that of a neutron, we replace the bare proton rms 
radius by an effective nucleon rms charge radius (ref. 6) which accounts 
for this difference. From reference 6, It Is 

- (N/Z) rl (5) 

where the bare proton rms radius Is rp ■ 0.87 fm (ref. 2) and the neutron 
rms radius Is r^ ■ 0.3359 fm (ref. 8). In equation (5)> N Is the neutron 
number and Z the proton number for the nucleus unde>' consideration. 
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RESULTS 


Abrasion cross sections for neon-carbon collisions, using equation 
(1), are listed In Table II as a function of nucleon charge radius, r^, 
and Fermi momentum, kp. The values of r^ correspond to the bare proton 
(0.87 fm) and the effective nucleon radius (0.806 fm) obtained from 
equation (5). The values for kp correspond to the Infinite matter value 
(1.36 fm“M» the value for the composite system (A ■ Ap + Ap where 
kp » 1.23 fni"M» «f»d a representative value, from reference 5, idilch 
Includes drift momentum effects (0.7 fm"^). In order to compare these 
results with the experimental data of Stevenson et a1 (ref. 3), It 1$ 
necessary to convert the abrasion cross sections Into relative 
probabilities for the formation of a particular residual projectile 
fragment mass, Ap. The results of this procedure, which Is described In 
detail In reference 2, are listed In Table III and displayed In figures 1 
and 2. Also displayed In the figures are the experimental data of 
reference 3. 

In figure 1 are displayed the results obtained for Ne+C In reference 2 
(rc • 0.87 fm, kp ■ 1.36 fin'M this work (r^ ■ 0.806 fm, 
kp ■ 0.7 fm*M compared with the experimental data (ref. 3). The 
agreement between experlmenc and the results of this work Is excellent. 

From Table III, analysis of the relative probabilities Indicates that there 
Is essentially no dependence of the relative probability on kp. Note 
that Table II shows oabs decreasing as kp decreases. All values of 
oabs listed In Table II, however, are In excellent agreement with the 
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experimental value of 1040±60 mb (ref. 9). From figure 2 we see that the 
Improved agreement with the experiment Is due to Incorporating the neutron 
charge distribution differences (from eq. (5)) Into the effective nucleon 
charge distribution. The Fermi momentum was kp > 0.7 fm~^ for both 
curves. 

Figure 3 displays results obtained for Ne + Mo with r^ " .00 fm and 
kp « 0.7 fm~^. In general, the agreement Is quite good except when 
Ap £ 4 where the theory overestimates the relative probabilities. The 
theoretical curve was determined using ^^Mo as a target. The experiment 
data were obtained for natural molyboenun which has 7 stable Isotopes 
(A « 92 , 94 , 95 , 96 , 97 , 98, 100) of roughly compar<*ble abundance (9.04 
percent to 23.78 percent). Unfortunately, experimental charge distribution 
data (ref. 8) are not available for all stable molybdenum Isotopes so that 
a more exact theoretical analysis Is not possible. 
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CONCLUDING REMARKS 


By utilizing an effective nucleon charge distribution, rather than the 
bare proton distribution, to account for differences in the charge 
distributions within the proton and neutron, improved agreement between .he 
predicted abrasion cross sections and recent experimental data were 
obtained. These findings also confirm the sensitivity of the abrasion 
results to the assumed nuclear distribution found in reference 2. In 
reference 2, the need for Pauli effects to be included were also clearly 
demonstrated. In this work, however, we find that once Pauli effects are 
included, the abrasion results are relatively insensitive (at 2.1 
GeV/nucleon) to the actual value of Fermi momentum used, as long as it is 
physically realistic (less than the infinite matter value). Further 
confirmation of these findings and additional Improvements to the theory 
will require additional experimental data. 
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Table II 


Abrasion Cross Sections for ^®He + X 

(Incident kinetic energy Is 2.1 GeV/nucleon). 


Number of Abraded 
nucleons, n 

Abrasion Cross Sections, mb 

rc « 0.87 fin 

rc«.806 fin 


kfl.36 fin-‘ 

kp-1.23 

kp-0.7 fin-^ 

kp-0.7 fin-‘ 

1 

248 

241 

241 

252 

2 

134 

131 

131 

136 

3 

96 

93 

93 

96 

4 

76 

74 

74 

77 

5 

64 

63 

63 

65 

6 

57 

56 

56 

57 

7 

52 

50 

50 

52 

8 

48 

47 

47 

48 

9 

45 

44 

44 

46 

10 

43 

42 

42 

44 

11 

42 

41 

41 

42 

12 

40 

39 

39 

40 

13 

37 

37 

36 

36 

14 

33 

33 

32 

31 

15 

26 

26 

25 

24 

16 

18 

18 

16 

16 

17 

10 

10 

9 

8 

18 

4 

4 

4 

3 

19 

1 

1 

1 

0.8 

20 

0.1 

0.1 

0.1 

0.1 

oabs 

TUT! 
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mr~ 

K57? 
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Table III 


Relative Probabilities for Fonaatlon of Projectile 
Fragnent Ap for Ne -t- C 


Relative Probability 


Af 


- 0.87 fm 




kp - 1.36 fii -‘ 

kp - 1.23 fii -‘ 


kp - 0.7 fi -‘ 

19 

.131 

.130 

.131 

.133 

18 

.141 

.140 

.141 

.143 

17 

.101 

.100 

.101 

.102 

16 

.080 

.080 

.080 

.081 

15 

.068 

.068 

.068 

.068 

14 

.060 

.060 

.060 

.060 

13 

.054 

.054 

.055 

.055 

12 

.050 

.050 

.051 

.051 

11 

.048 

.048 

.048 

.048 

10 

.046 

.046 

.046 

, V . 

9 

.044 

.044 

.044 

. 04 s 

8 

.042 

.042 

.042 

.042 

7 

.039 

.040 

.039 

.038 

6 

.035 

.035 

.034 

.033 

5 

.027 

.028 

.027 

.025 

4 

.019 

.019 

.019 

.016 

3 

.010 

.011 

.010 

.009 

2 

.004 

.004 

.004 

.003 

1 

.001 

.001 

.001 

.0008 

0 

.0001 

.0001 

.0001 

.0001 
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Figure 1. Abrasion results for neon projectiles colliding 
with carbon targets, as predicted by this work and the 
previous abrasion model (ref. 2), compared with experiment. 
Incident kinetic energy is 2.1 GeV/nucleon. 
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Figure 3. Abrasion results for neon projectiles 
colliding with molybdenum targets. Incident kinetic 
energy is 2.1 6eV/nucleon. 


